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We provide a constructive algorithm to find the best sepa- 
rable approximation to an arbitrary density matrix of a com- 
posite quantum system of finite dimensions. The method 
leads to a condition of separability and to a measure of en- 
tanglement. 
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Entanglement and nonlocality are some of the most 
emblematic concepts embodied in quantum mechanics 
0. The non-local character of an entangled system is 
usually manifested in quantum correlations between sub- 
systems that have interacted in the past but are not 
longer interacting. Furthermore, these concepts play a 
crucial role in quantum information theory . 

From a formal point of view, a state of a composite 
quantum system is called "inseparable" (or "entangled" ) 
if it cannot be represented as a tensor product of states 
of its subsystems. On the contrary, a density matrix p 
describes a separable-state if it can be expressed as a finite 
H sum of tensor products of its subsystems: 

p. = ^2pi(Pi®Pi f ~®p?y, i>k>o (i) 

i 

where pf,pf,.., pf are density matrices describing sub- 
systems A,B, ...,N, respectively and YliPi = 1- Thus, 
separable states are those than can be produced by N dis- 
tant observers (Alice, Bob,..,Norberto) that prepare their 
states (pf , pf , pf) independently, following common 
instructions (pi) from a source [Q. Let us, for the mo- 
ment, restrict ourselves to binary composite systems, i.e. 
TL = TLa ®7~Lb- Using the spectral decompositions of pf 
and pf it is easy to rewrite the Eq. (|l|) in the form 

p s =J2 X « P <* l>A a >0; 5Z A « = 1 ' ( 2 ) 

a a 

where a is a multi-index running over all distinct eigen- 
vectors of the matrices pf ® pf, and P a are projec- 
tors onto product states, i.e. P a = |e, /)(e,/| (where 
|e) e Ha and |/) G Hb)- Separable states, p s , are thus 
mixtures of product states and as such their correlations 
are purely classical. 

The distinction between entangled and separable states 
is well established for pure stated: entangled pure states 
do always violate Bell inequalities || . For mixed states, 
however, the statistical properties of the mixture can hide 
the quantum correlations embodied in the system, mak- 
ing thus the distinction between separable and entangled 



enormously difficult |^^] . Besides the importance of the 
subject from a fundamental point of view, this distinc- 
tion has also important consequences for quantum infor- 
mation theory. Consider, for instance, Werner's family 
of entangled mixed states j^] , that does not violate any 
kind of Bell inequalities but, nevertheless, can be used 
for quantum teleportation 

Recently, a first step in such distinction has been done 
by Peres |J] and the Horodecki family [p|JIo|. They have 
formulated two necessary conditions to characterize sepa- 
rable density matrices. The first condition Q| states that 
if a matrix p is separable, then its partial transposition 
(with respect to subsystem A, or B) must be a density 
matrix, i.e. must have non-negative eigenvalues: 

P = Ps ^P Tb = (p Ta T >0. (3) 

This can be easily grasped from the representation (||) of 
separable matrices, since the partial transposition with 
respect to system B, amounts to replacing P a by P^ B = 
|e, /*)(e, /*|, so that evidently 

P TB =E A «l e '/*>< e '^°- ( 4 ) 

a 

This condition is sufficient to guarantee separability oly 
for composite systems of dimension 2 x 2 or 2 x 3. 

The second necessary condition || states that if p = 
p s , then there exist a set of product vectors V = {\ei,fi}} 
that spans IZ(p) and at the same time V T2 — {\ei, /*)} 
spans lZ(p T2 ) where IZ(p) denotes the range of p, i.e. 
the set of all \ip) G Ti, for which 3 \ <j>) G H such that 
= p\4>)- From the representations (||) and (Q) we 
see that if a set of product vectors {\ei, fi)} spans IZ(p), 
it immediately follows that the set of product vectors 
{\ei, /*)} also spans 1Z(p TB ). In general, both conditions 
are not equivalent. In particular, when the dimension 
of IZ(p) is equal to the dimension of lZ(p T2 ), the second 
condition may not be sufficient to ensure separability. 

Finally, let us point out, that for a density matrix 
which is known to be separable, only if dim[7i] < 6 there 
exist an algorithm for decomposing it according to Eq. 

In this Letter we address this last point and provide a 
constructive way of finding such an algorithm regardless 
the (finite) dimension of the composite system. That im- 
mediately leads to a necessary condition for separability. 
Furthermore, we shall demonstrate that any inseparable 
mixed state in C 2 <£> C 2 can be decomposed in a separable 
matrix and just a single pure entangled state, providing 
thus a novel characterization of the "entanglement" of 
any inseparable state. 
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The idea behind the algorithm relies on the fact that 
the set of separable states is compact. Therefore, for 
any density matrix p there exist a "maximal" separable 
matrix p* which can be subtracted from p maintaining 
the positivity of the difference, p — p* > 0. Let us express 
the above idea in a more rigorous way: 

Theorem 1 For any density matrix p (separable, or not) 
and for any set V of product vectors belonging to the 
range of p, i.e |e, /) £ T^(p) there exist a separable (in 
general not normalized) matrix 



(5) 



with all A Q > 0, such that 5p = p — p* > 0, and that p* 
provides the best separable approximation (BSA) to p in 
the sense that the trace Tr(<5p) is minimal (or, equiva- 
lently, Trp* < 1 is maximal). 

The proof of the theorem is simple, and the whole art 
is, of course to construct pt. Let us consider all separable 
matrices p s of the form (|5|) that we can subtract from p 
maintaining the non-negativity of the difference 5 p. Ob- 
viously, the trace of p s must be smaller than one, since 
< Tr(5p) = 1 — Trp s . The set of such matrices is de- 
termined by the set of possible A a > for which 5p > 0, 
and < Trp s — J2 a A a < 1. This set is closed (in any 
reasonable topology) . The set of all possible traces of p s 
is bounded from above, so it must have an upper bound, 
say 1 — e; ergo because of the compactness of the set of 
all p s , there exist a matrix p* in this set with the maxi- 
mal trace, equal to 1 — e. That implies that although the 
matrix p* s \V] depends on the choice of the set V, and by 
expanding V we can construct better separable approxi- 
mations to p (i.e. for V D V, Trp*[V] > Trp*[V]), it is 
generally sufficient to take V C S large enough to obtain 
already the maximal possible trace Trp*[F] = Trp*[S] 
(where S is the set of all \e,f) S TZ{p)). The latter 
statement indicates also that although typically the BSA 
matrix Ps[V^] is not unique, its trace is. Nevertheless, 
for C 2 <g) C 2 composite systems we shall demonstrate that 
Ps[V] is also unique. 

As an obvious consequence of Theorem 1, we obtain a 
necessary and sufficient condition for separability: 

Condition 3 A density matrix p is separable iff there 
exist a set of product vectors V C TZ(p), for which the 
best separable approximation to p, /0*[V] has the trace 1. 

The proof is again simple: The necessity of the cond3 
follows directly from (||). From the fact that Sp = p — 
P* s > 0, and Tr<5p = 1 — 1 = Q, we obtain Sp = 0, or 
equivalently p = p*. 

Before we discuss the procedure of construction of the 
matrix p*, let us to introduce two concepts which shall 
play a crucial role in what it follows. 

Definition 1 A non-negative parameter A is called max- 
imal with respect to a (not necessarily normalised) den- 
sity matrix p, and the projection operator P — \ip){ij}\ iff 



p — AP > 0, and for every e > 0, the matrix p — (A + e)P 
is not positive definite. 

The maximal A determines thus the maximal contribu- 
tion of P that can be subtracted from p maintaining the 
non-negativity of the difference. In the following we will 
apply the above definition to projections onto product 
vectors, i.e. \tp) — |e, /). The following lemma charac- 
terizes a single maximal A completely: 

Lemma 1 A is maximal with respect to p and P = \ip){ip\ 
iff: (a) if \tp) K(p) then A = 0, and (b) if |</>) € K{p) 
then 



< A = 



f 



(6) 



Note that in the case (b) the expression on RHS of Eq. 
@ makes sense, since € H(p), and therefore there 
exists € H(p) such that = p\^)- Let us observe, 
that for any \(f>) the Schwartz inequality implies that 

{<t>\p\<t>) = m^p^m 2 < (4\p\mfyi>)- (7) 

That proves that for every \<j>), (4>\p— (ip\l / plip) -1 P\</>) > 
0, i.e. p — AP > 0. Since on the other hand, (p — 
AP))*) = for |*> = ±\ip), thus for every e > 0, 
(A + e)P]m = -eA~ 2 < 0. This proves that A given by 
expression (||) is indeed maximal. 

Definition 2 A pair of non- negative (Ai,A2) is called 
maximal with respect to p and a pair of projection opera- 
torsP! - P2 = \ih){ih\ iff p-A 1 P l ~A 2 P 2 > 0, 

Ai is maximal with respect to p — A2P2 and to the pro- 
jector Pi, A2 is maximal with respect to p — A\Pi and to 
the projector P2, and the sum Ai + A2 is maximal. 

The maximal pair (Ai, A2) determines thus the maxi- 
mal contribution of A1P1 + A2P2 that can be subtracted 
from p maintaining the non-negativity of the difference, 
and that has a maximal trace, Tr(AiPi + A2P2) = 
Ai + A 2 . 

Lemma 2 A pair (Ai,A2) is maximal with respect to 
p and a pair of projectors (Pi,P2) iff: (a) if \1fi1), \ip2) 
do not belong to IZ(p) then Ai = A2 = 0; (b) if |V , i) 
does not belong to TZ(p), while ^2) S TZ(p) then Ai = 
0, A 2 = (c) if IVi), 1^2) € Tl{p) and 

(Viil/p|^2> = then Ai = (i/ji\l/p\^i), * = 1,2; (d) 
finally, if |V>i>, ^2) € K(p) and (^1/ p\ifj 2 ) f then 



Ai - «Vs|l/p|V2> 
A 2 = ((^i|l/p|Vi) 



K^ll/pl^D/P, (8a) 

|^l|l/p|^2>|)/P, (8b) 



where D = / p^) {^ 2 \1 / p\^ 2 ) - |<^i|l//#2>| 9 . 

The proof of (a) and (b) is the same as the proof of 
Lemma 1. In the case (c) observe that (p— AiPi)^ 1 \tp 2 ) = 
p" 1 !^), (p- A 2 P 2 y 1 \-ipi) = 1 , so that maximal- 
ly of Ai implies automatically that Aj = (^ilp -1 !^), i = 
1,2. Finally, in the case (d) we get (p — A 2 P2) _1 IV'i) = 
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p" 1 ^} +Bp- 1 \4> 2 ), with B = A 2 (^ 2 |l/p|V>i)/D. The 
maximality of Ai assures then automatically the maxi- 
mality of A 2 provided 

1 - A 1 (i; 1 \l/p\<ip 1 ) - A 2 <V> 2 |l/p|^ 2 ) + AxA 2j D - 0. (9) 

Maximizing the sum Aj + A 2 with the constraint (^), we 

arrive after elementary algebra at Eqs. (||). 

We can now formulate the basic theorem of this paper: 

Theorem 2 Given the set V of product vectors |e, /) G 
IZ(p), the matrix p* = ^2 a A a P a is the best separable 
approximation (BSA) to p iff a) all A Q are maximal with 
respect to p a — p — ^2 a i^ a A a iP a t, and to the projec- 
tor P a ; b) all pairs (A a , Ap) are maximal with respect to 
Pap = p - J2 a '^ a ,p A Q <P Q ', and to the projection opera- 
tors (P a ,Pp). 

Let us prove now that maximizing all the pairs 
(A a ,Ap) with respect to p a p = p- J2 a '^ a ,p ^-a'Pa' , 
(P a ,Pp) is a necessary and sufficient condition to sub- 
tract the "maximal" separable matrix p* from p. Obvi- 
ously, if p* is the BSA then all A Q , as well as all pairs 
(A Q , A^) must be maximal, since otherwise maximalizing 
A Q , or the sum A Q + A^ would increase the trace of p*, 
maintaining non- negativity of p — p* . 

To prove the inverse, assume that the total number 
of a's is K, and that p* has all pairs of A's maximal. 
Consider matrices p s — ^ a X a P a in the vicinity of p*, 
for which all individual X a are maximal, i.e. p s belong to 
the boundary of the set Z of all separable matrices such 
that p — p s > 0; A Q 's lie thus on a (K — l)-dimensional 
manifold, defined through a constraint, 



f(X 1 ,...,X K ) = 0. 



(10) 



Maximality of (A Q , A^) implies that (X a + Xp) has a max- 
imum at X a .p = A a .p under the constraint JlOj), and for 
all 7 a, f3] A 7 = A 7 which implies (df /dX a \x=A) = 
{d f J dXp\\—fC) ■ Using this identity for sufficient number 
of pairs we get that (<9//i9A q ,|a=a) = const for all a. 
That is equivalent to the fact that the gradient of Tr(p s ) 
under the constraint (|l^) vanishes for p s = p* . The trace 
of p s has thus either a local maximum, or a minimum, 
or a saddle point at A = A. The two latter possibilities 
cannot occur, since the trace is maximal with respect to 
all pairs of A's, and since the set Z is convex ( i.e. if 
Ps,p' s G Z then ep s + (1 - e)p' s G Z for every < e < 1). 
For the same reason of convexity, the local maximum at 
p* must be a global one, i.e. there cannot exist two ma- 
trices p*, and p*, which both provide local maxima of 
the trace, and have Trp* 7^ Trp*; ergo p* is the BSA, and 
any other matrix p* which has all pairs of A's maximal, 
must have the same trace as p*. 

In any case, we have shown that any density matrix 
p of composite system Ji can be decomposed according 
to p = p* + Sp where p* is a separable matrix (in gen- 
eral not normalized) with maximal trace. Let us analyze 



such decomposition in more detail. All the information 
concerning "inseparability" is included in the matrix 5 p. 
If it does not vanish, i.e. if p is not separable, its range 
TZ(Sp) cannot contain any product vector. We have ob- 
served that, quite typically, if Sp is a sum of projections 
onto a set of linearly independent entangled states, then 
there exist product vectors that belong to lZ(Sp), whose 
contributions can be single out increasing Trp s *. The 
reason is that, for instance, the set of all product vec- 
tors in the Hilbert space H of dimension N x M spans 
a (N + M — l)-dimensional manifold, which generically 
has a non- vanishing intersection with linear subspaces of 
H of dimension equal or larger than (N — 1) X (M — 1). 
The above statement implies that for N = M = 2, Sp is 
a simple projector onto an entangled state. 

As an immediate consequence, we obtain that any den- 
sity matrix p in C 2 ® C 2 has a unique decomposition in 
the form: 



p = Ap s + (l-A)P e ; AG [0,1] 



(11) 



where p s is a separable density matrix (normalized), 
P e denotes a single pure entangled projector [P e = 
|4 , e )(5 , e |), and A is maximal. Any other decomposition 
of the form p = Ap s + (1 — A)P e with A G [0, 1] such 
that p s ^ p s , necessarily implies that A < A. If not, that 
is, if A = A for p s ^ p s , it follows from Ref. jll| that 
for P e 7^ P e , we can always find projectors onto product 
states in the plane formed by P e and P e and therefore in- 
crease A, which is impossible since A is already maximal. 

The decomposition given by expression (|ll]) leads 
straightforwardly to an unambiguous measure of the en- 
tanglement for any mixed state p (in C 2 (g> C 2 ): 



£(p) = (l-A)x£(|* e » 



(12) 



where E(\^ e )) is the entanglement of its pure state ex- 
pressed in terms of the von Neumann entropy of the re- 
duced density matrix of either of its subsystems Jl^] : 

£(|* e » = -Trp A log 2 p A = -Trp B log 2 p B (13) 

where P{a,b} = ^{B.AjP- This measure of entanglement 
is clearly independent of any purification or formation 
procedure ]l2] , p^[ . 

Let us illustrate with an example the ideas stressed 
in the paper. Consider a pair of spin-^ particles in an 
impure state consisting of a fraction x of the singlet and a 
mixture in equal proportions of the singlet and the triplet 
||. This state is described, in the computational basis, 
by the density matrix: 



Pw(x) 



/ 1 — x 
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< x < 1 (14) 



For this case Eq.(||) is sufficient to ensure separability: 
p w is separable if x < 1/3 and inseparable otherwise. 
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Nevertheless, we use our procedure to check the separa- 
bility and to obtain the decomposition of p given by Eq. 
( |i"l| ) for different values of x. 
For each given set V, we first construct the matrix 



Pirn 



v 



(15) 



with the A' maximized pairwise, according to the defi- 
nitions fl4|| . When the numerical convergence has been 
achieved we obtain Sp = (p w — /o*[V]) and compute its 
trace. Typically, we observe that: (a) only very few pro- 
jectors P a of each set V contribute to the matrix /0*[V], 
and (b) if the set V is large enough (i.e > 300), the results 
become independent of the chosen set. 

The results are presented in Fig.l, for a set of 100, 
200 and 500 P a -projectors randomly chosen. Each point 
represents the corresponding value of Tr (Sp) for a given 
p w (x). The vertical line indicates the condition of sepa- 
rability, derived from Eq.(|). For x < 1/3 , Tr (Sp) = 
indicating that p w is separable. At x ~ 1/3, a clear 
"phase-transition" occurs, and the value Tr(Sp) ^ 0, in- 
dicating thus the non-separable character of the state. 
Therefore, our numerical results reproduce accurately the 
conditions of separability derived from Eq.(0). 



x > 1/3. A measure of the entanglement of p w is, there- 
fore, naturally provided by the value of the corresponding 
A, i.e. E(p w (x)) = (1 — A(a;))-ebits, since the singlet has 
a value of entanglement of 1 e-bit (see (Eq.(|ilj|)). This 
measure does not coincide with other measures of the 
entanglement of p w (lj,[ll|. A further analysis of this 
entanglement measure will be presented elsewhere. 

Summarizing, we have presented a method to construct 
the best separable approximation to an arbitrary den- 
sity matrix of a composite quantum system (of arbitrary 
dimensions). The method provides a necessary condi- 
tion for separability of a density matrix. Furthermore, 
for composite systems of dim[W]=4, it also provides with 
unambiguous measure of the entanglement of its non sep- 
arable states. 
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FIG. 1. The best separable approximation to a Werner 
state p m . We plot the the value Tr(<5p) for the matrix p w 
characterized by the fraccion of the singlet x (see(|14)). The 
vertical line indicates the separability border (Eq.pl)). (The 
numerical precision of the algorithm is set to 10~ , so that 
Tr(6p) must be > 10" 4 ). 

Let us now analyze the "inseparability" properties of 
p w . The matrix Sp when it does not vanish, i.e. for 
x > 1/3 corresponds to the projector onto the maximally 
entangled singlet = l/v%(| T!) - I IT))- Thus, a 

Werner state of the type p w can always be decomposed 
as: 

p w (x) = X(x) Ps + (1 - A(a;))|¥-)<*-| (16) 
with A = 1 for a; < 1/3 ( p w = p s ), and < A < 1 for 
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